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Abstract 

A dilute distribution of magnetic impurities is assumed to be present in doped graphene. We 
calculate the interaction energy between two magnetic impurities which are coupled via the indirect- 
exchange or Ruderman-Kittel-Kasuva-Yosida (RKKY) interaction by the doped conduction elec- 
trons. The current model is a half-filled AB-laitice structure. Our calculations are based on the 
retarded lattice Green's function formalism in momentum-energy space which is employed in linear 
response theory to determine the magnetic susceptibility in coordinate space. Analytic results are 
obtained for gapped graphene when the magnetic impurities are placed on the A and B sublattice 
sites of the structure. This interaction, which is important in determining spin ordering, has been 
found to be significantly different for AA and BB exchange energies in doped graphene due to the 
existence of an energy gap, and is attributed to a consequence of the local fields not being equal on 
the A and B sublattices. For doped graphene, the oscillations of all three RKKY interactions from 
ferromagnetic to antiferromagnetic with increasing Fermi energy is significantly modified by the 
energy gap both in magnitude and phase. Additionally, the AB exchange energy may be modified 
by the presence of a gap for undoped graphene but not for doped graphene due to the dominance 
of doped conduction electrons. 



PACS numbers: 75.30.Hx;75.10.Lp 



I. INTRODUCTION 



Recently, there has been some attention given to the effects arising from the Ruderman- 
Kittel-Kasuva-Yosida (RKKY)^ or indirect-exchange interaction between spins via the 
host conduction electrons of monolayer free standing graphene, a two-dimensional honey- 
comb network of carbon atoms^^. Such local moments may arise near extended defects. 
Interest in the RKKY interaction is due partially to the fact that it determines spin ordering. 
The indirect exchange interaction between local magnetic moments is generally determined 
by electron excitations near the Fermi level. However, although the RKKY interaction has 
been considered in intrinsic graphene by several authors, little attention has been given to its 
role when a gap 2E g is opened at the two inequivalent K and K' points in momentum space 
between the valence and conduction bands^HEI. This effective band gap may be generated 
by spin-orbit interaction^-. It may also arise when monolayer graphene is placed on a sub- 
strate such as ceramic silicon carbide^ or graphite 14 or dynamically when it is irradiated by 
circularly polarized light ^H^. Depending on the nature of the substrate on which graphene 
is placed or the intensity and amplitude of the light, the gap may be a few meV or as large 
as an eV. In general, the energy gap is attributed to a breakdown in symmetry between the 
sublattices caused by external perturbing fields from the substrate or photons coupled to the 
A and B atoms. Furthermore, there still exists a small band gap 2E so ~ 1CT 3 meV due to 
spin-orbit coupling! 18 * 19 * even in the absence of a substrate or an external laser field. Despite 
the formation of an energy band gap, corresponding to a metal-insulator transition, for the 
half-filled bipartite lattice, the interaction between atoms on the same lattice has still been 
suggested to be ferromagnetic and antiferromagnetic between atoms on different sublattices 
in undoped graphene^, even though the energy dispersion is no longer linear and isotropic 
as it is for gapless Dirac electrons. For undoped graphene, we show that the existence of 
an energy gap can significantly modify the magnitude of the magnetic interaction between 
impurities on the sublattices. More interestingly, for doped graphene, we demonstrate that 
the energy gap can drastically affect the nature of such a magnetic interaction. 

In Ref. [10] , with a small gap at the Fermi level due to spin-orbit interaction, the graphene 
layer was assumed to be undoped. The present paper explores quantitatively the role played 
by massive relativistic Dirac particles in both doped and undoped graphene on the RKKY 
indirect-exchange energy. In other words, how the magnetic effects depend on the chirality of 
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the electron states for doped graphene with an energy gap. Analytic results are obtained for 
gapped graphene when both impurities are located on sites belonging to the A sublattice, 
both on the B sublattice and when one impurity is on an A while the other is on a B 
sublattice site. We demonstrate numerically the large asymmetry between the AA and BB- 
exchange-interaction energy for gapped graphene. Our closed-form analytic results make it 
convenient for our predictions to be compared with experimental data once they become 
available. 



The remainder of this paper is organized as follows. In Sec. II we obtain the retarded 



Green's function matrix elements on the A and B sublattices for gapped graphene. In 



Sec.|III[ we use linear response theory to calculate the magnetic susceptibility by making use 
of our derived results for the Green's functions on the bipartite sublattices. The interaction 
energy of two spins at lattice sites and r' u , where /i, v — A, B, is then deduced for both 



undoped and doped gapped graphene. We conclude our paper in Sec. IV with some remarks 



II. GREEN-FUNCTION FORMALISM FOR GRAPHENE 



The Hamiltonian for monolayer gaped graphene in the absence of impurities has the form: 



W(k) 



where 



E n 



-7 /i (k) -E g 



(1) 



h (k) = exp (ik y a) + exp 



i\/3ak x iak v 



+ exp 



i\/3ak x iaky 



(2) 



with the carbon-carbon distance a = 1.14 A, 7 = —3.0 eV and k = (k x , k y ). Here, 2E g is the 
energy gap generated by some means, possibly by a substrate or circularly polarized light, as 
discussed in the Introduction. The retarded Green's function matrix in wave vector-energy 
(k, E) space is given by 
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with 



G AA (k,E) = ^-^( Eg + E) , 
G BB (k,E) = ^-^(- Eg + E) , 
Gab (k, E) = G BA (k, E) = ^l_ 7o />o(k) , 

(4) 



where we introduced the energy dispersion equation E^ = 7 q hoik) h^ik) + E 2 g . 

We may now transform the Green's function to the real space representation with use of 



G(v,r',E) = J d 2 k G(k, E) exp [ik ■ (r - r')] 



S 

4^2 



1 st BZ 

6 

£ / d^dfc,, G(Ki + k, E) exp [i (K, + k) • (r - r')] , (5) 

i=i "* 

where the summation is carried out over the six corners of the Brillouin zone (BZ), so 
that we may make the approximation E^_. +k = (7/c) 2 + E 2 with 7 = 3a7 /2. In that 
approximation, we also have 70/10 = jk. Additionally, in our notation, the area of the unit 
cell is S = 3V3a 2 /2. 

Using cylindrical coordinates, the matrix elements of the retarded Green's function may 
be explicitly written as 

2tt k c 

G AA (r A , r' A ; E) = ^ cos^ • (r A - r' A )} J d9 J dk kG AA (k, E) e^A-r'j , 



2-7T k c 

G BB (r B ,r' B ;E) = cos^ • (r B - r' B )\ J d6 J dk kG BB (k, E) ^b-V b ) ? 



G AB {v A ,v' B -E) = G BA (r' B , r A ; E) 

= ^p{ e ' Kl(rA ~ r ' s) J d0 J dk k(-i sin 9 - cos e)G AB (k,E)e l ^ rA - r '^ 





2-7T k c 



+ e-*™fmf*«-i*>e + m WM**^, (6) 
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where k c is a cut-off wave vector whose existence is due to the fact that the energy band 
structure used in the calculation is only valid for a limited range of the first BZ near the K 
points. However, it is also allowed in such an approximation to take the limit k c — > oo from 
a mathematical point of view. In the above expression, we have Ki = (47r/3-\/3a, 0). 

We now make use of the following identity involving the series of the first kind of Bessel 
function J n (x) 

oo 

e *k. (l> -0 = Mklr ^ _ + 2 i n J n {k\v, - r'J) cos^,,,,,) , (7) 

n=l 

where (fx, v) = (A, B) and (p rfl ,r' u is the angle between k and r M — r' v . We may calculate the 
above integrals and simplify the real-space Green's function matrix elements as 



G AA {v A ,v' A ;E) = cos^ ■ {r A - v'j^F^v A -v' A \,E) , 
G BB {v B y B -E) = cos[K 1 .(r B -r' B )]F 2 (|r B -r' B |, J E?) , 

G AB (r A ,r' B ;E) = G BA {v' B , r A ; E) = sinfKx • (r A - r' B ) + a rA y g } F 3 (\r A - r' B |, E) . (8) 

In this notation, cos{ip r r iJ) = cos(# — a r r '„), where 6 is the angle which k makes with the 
x-axis, and a r y v is the angle between — r' u and the x-axis. Three position-dependent 
complex functions Fj(|r M — r'J, E) for j = 1, 2, 3 are defined as follows: 
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where K m (x) for integer m is the modified Bessel function. In addition, we have used the 
formula for the m th order Hankel transform 



oo 



/x m 
xdx -. J m (xE) = a m K m (aE) , (10) 
x z + cr 

o 

where a is an arbitrary constant. 



III. INDIRECT-EXCHANGE INTERACTION ENERGY 

Let us now assume that we have two magnetic impurities with spins Si and S 2 ) located at 
r M and r' u , respectively. According to linear response theory, the energy needed to exchange 
(mediated by the Dirac electrons) their positions may be written in the matrix form 

V^d^-i/J) = ^(1^-^1)81.82, (11) 

where the exchange-integral matrix is proportional to the matrix of spin-independent sus- 
ceptibility 

J»*(\r»-<\) = A ^X,»(\r»-rl\) • (12) 

Here, A is the contact interaction between the magnetic impurities and the Dirac electrons. 
The magnetic susceptibility may be expressed in terms of the Green's functions in the 
standard RKKY form^. 



A. Undoped Graphene with a Gap 

For undoped graphene, we take the Fermi energy Ep = 0. For this case, we obtain the 
matrix of spin-independent susceptibility 

-Eg 

X^lr, - r'„|) = -1 J dE Im [G(r M , r' u ; E) G(r' u , r^; E)] . (13) 



-oo 



Making use of the expressions for the Green's functions from the preceding section, it is a 
simple matter to obtain the following matrix elements: 
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^9 

X aa{\va-v' a \) = -^(l + cospKx-^-i^)]) Im j dE [F^ - r' A \, E)} 2 , 

— oo 
-Eg 

X bb{\vb-v' b \) = ~(1 + cos[2K 1 -(t b -t' b )]) Im J dE [F 2 (\r B - r' B \, E)} 2 , 
Xab(\ta - r' B \) = -- (cos[2K! • (r A - r' B ) + 2a rA y B ] - l) 



x 



^9 

Im j dE[F z {\v A -v' B \,E)] 2 . (14) 



We first consider the case when both impurities are located on A atomic sites. That is, we 
need to calculate 



V^.Ru) = Im J dE[F l (\r A -r> A \,E)Y 



^JdE(E + E,f J h^5|, - r'A iVo | r , _ r ; 4 | 



2 

37o 



J dE (E + Eg) 2 J ( ^ E ° 

— oo 
-Eg 

J dE (E + E g ) 2 Jo (~\^E* - El R AA ^j N E 2 ~ E] R AA ^j , (15) 



where N (x) is the Neumann function. For convenience, in the last expression, we measure 
energy in units of 70 and distance in units of R AA = \r A — r' A \/a. In the gapless case, E g = 0, 
we are able to calculate the integral analytically, obtaining 



u 

Vl{R AA ) = V 1 (0,R AA ) = Im J dE [F^ta - r' A \,E)f 

— oo 

00 
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, dx x 2 J (x) N (x) = —5 — . (16) 

A l0 R\ A J oy } oy } 6A lo R AA y } 







The last integral may be evaluated through 
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lim dxx exp(— xqx) Jq (x) Nq (x) 

x ^0 J 

o 

lim 4 (4 + 3X ° } £ (- X »/^ ~ 2 ( 8 + 6x o + X § ^(-^o/ 4 ) _ 1 ( 17) 
where /C(x), are the complete elliptic integrals of the first and second kind, respectively. 

In a similar fashion, we have for the AB exchange interaction 



V 3 (E g ,R AB ) = Im J dE[F 3 (\r A -r' B \,E)] 2 = -^- J dE (^J E 2 - 

— oo — oo 

x J x \-\^E* - E*Rab\ JVi \-\^E 2 - E*R AB \ , (18) 

where Ni(x) is also the Neumann function. For E g = 0, the integral may be evaluated 
analytically, yielding 



lim dx x exp(— xqx) J\ (x) Ni (x) 

x ^0 ' 







-4 (-4 + xg) g(-xg/4) - 4 (4 + xg) K(-sg/4) = 



so that we have 



u 

V^R AB )=V 3 (0,R AB ) = Im j dE[F 3 (\r A -r' B \,E)] 2 = - 6A 2 J R3 . (20) 

— oo 

The effects of doping on the indirect-exchange interaction for a graphene layer without 
a gap have been studied in detail. In this case, a sign change of the indirect-exchange 
interaction was discovered when the doping concentration was varied. There is a switching 
from inverse cubic to inverse square power law for the spatial dependence.^ However, the 
doping effect for a graphene layer with a gap remains largely unexplored. To demonstrate 
the doping effect for gaped graphene, we write Vj(E g , R^ u ) = Vj(R^ u ) + AVj(E g , R pv ) with 
[A, v — A, B and j = 1, 2, 3. 
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In Fig.[TJ we present our calculated results for Vj(E g , R^) (left panel) and AVj(E g , R^) 
(right panel) in an undoped graphene layer as functions E g with Raa = Rbb = Rab = 0.5. 
From the left panel of this figure, we find that the positive exchange interactions V\ and 
V2 between two intra-sublattice impurities (intra-SIs) decrease with E g while the negative 
exchange interaction P3 between two inter-sublattice impurities (inter-SIs) increases with E g 
at the same time. We further observe that the degenerate V\ and V2 at E g = 0, due to the 
presence of symmetry with respect to two sublattices, is now lifted for E g > 0. Compared 
with the results when E g = 0. We see from the right panel of this figure that the gap- induced 
changes AVi and AV2 become more and more negative with increasing E g . Additionally, 
the decrease of AVi is more rapid than for A7V This behavior reflects a modulation in 
the exchange interaction by an energy gap in graphene. Also, it suggests a possible gap- 
induced ferromagnetic-to-antiferromagnetic phase transition due to the intra-SI exchange 
interaction. In contrast, a distinctly opposite behavior is found in the inter-SI exchange 
interaction AV3. 

We display in Fig.[2] the calculated results for R 3 ^ u Vj(E g , i? M „) (left column) and 
R^ u AVj(E g , R^y) (right column) as functions of the distance R^ v between two magnetic 
impurities for various values of E g in undoped graphene. From the left column of this figure, 
we observe clear deviations of the exchange interactions Vi, V2 and V3 from 1/R? depen- 
dence (at E g = 0) with increasing E g (from curve-1 to curve-3). These results also indicate 
a possible gap-induced phase transition which may be driven by either the intra-SI or the 
inter-SI exchange interaction at a relatively large distance between two magnetic impurities 
in a graphene layer. Additionally, we find from the right column of this figure that changes 
in AVi and AV2 in the intra-SI exchange interactions drop more rapidly as a function of 
R^n when E g is increased. Again, a completely opposite behavior is observed in this case 
for AV3 which is associated with the inter-SI exchange interaction. 



B. Doped Graphene with a Gap 

We now turn to the case of doped graphene which has a finite Fermi energy Ep > 0. It 
is convenient to introduce a dimensionless f the Fermi energy as Xp = ^EpRaa/Io- The 



Fermi energy makes the integration limits in Eqs. (16) and (20) raised to include the part 
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J dE. In the case when E g = those integrals assume the analytical forms: 
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Here, (?™' n is the Meijer's G-function. It has oscillating behavior^ thus switching the nature 
of Jaa and Jab between ferromagnetic and antiferromagnetic. In the long distance limit 
Raa ~ Rab ~ Rbb-i we can use the mean field approximation as in the work of Sherafati 
and Satpathy^. 

In Fig.[3j we present our calculated results for Vj(E g , R^) (left panel) and AVj(E g , R^) 
(right panel) as well as functions of E g with Raa = Rbb = Rab = 1 and Ep — E g = 1 in 
doped graphene. The left panel of this figure shows that the intra-SI exchange interaction 
V\ (V2) decreases (increases) with E g while the inter-SI exchange interaction V3 is almost 
independent of E g at the same time. The anomalous increasing P 2 with E g is in strong 
contrast with the results shown in the left panel of Fig. [T] and is indicative of the band-filling 
effect on the gap-modulation of the indirect exchange interaction between two magnetic 
impurities. Moreover, from the right panel of this figure, we observe that the gap-induced 
changes AV\ and AV2 acquire a linear dependence on E g in the presence of doping, which 
is also in contrast with the results presented in the right panel of Fig. [T] where a nonlinear 
dependence on E g is obtained. 

Figure [i] presents the calculated results for Vj(E g , R^ v ) (left column) and AVj(E g , R^) 
(right column) as functions of Ef — E g with Raa = Rbb = Rab = 1 and various values 
of Eg in doped graphene. From the left column of this figure, as predicted before^ we find 
that the band-filling introduces an oscillation between positive (ferromagnetic) and negative 
(antiferromagnetic) signs in V\, V2 and TV In addition, the zeros (non-magnetic) in TV 
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Vi and V3 are shifted leftward by increasing the values of E g . From AVi, AV2 and AV% 
displayed in the right column of this figure, we further discover that the magnitude of such 
oscillations is enhanced significantly in the presence of an energy gap for a graphene layer. 

We exhibit in Fig.[5] our numerical results for R^ lll Vj(E g , R pu ) (left column) and 
R 2 v AVj{Eg, R pv ) (right column) as functions of the distance R^ u between two magnetic 
impurities with Ep — E g — 1 and various values of E g in doped graphene. From the left 
column of this figure, we find that both positive intra-SI exchange interactions V\ and V2 
decrease with R^ but the increase of E g suppresses (enhances) V\ (V2), respectively. The 
mixed 1/-R5L and dependence in VipiEg, R^), as predicted in Ref . [H] for doped 

graphene, is modified dramatically. This clearly demonstrates the gap modulation of the 
doping effect on the intra-SI indirect exchange interaction between magnetic impurities in 
graphene. Moreover, P3 is found to increase with Rab with a negligible effect from E g , 
which is accompanied by a switching from antiferromagnetic (negative) to ferromagnetic 
(positive) beyond a threshold distance between a pair of magnetic impurities in the inter-SI 
exchange interaction. From the right-hand column of this figure, we find that these exists a 
negative minimum (positive maximum) for the gap-induced change AVi (AV2) as a function 
of RjjLfjD respectively. Interestingly, the gap-induced change AV3 for the inter-SI exchange 
interaction shows up a sign change with increasing Rab for all chosen finite values of E g . 

IV. CONCLUDING REMARKS 

The RKKY interaction between a pair of magnetic impurities located on lattice sites 
of monolayer half-filled gapped graphene is theoretically investigated based on the lattice 
Green's function formalism. In contrast with the case for gapless monolayer graphene, 
the RKKY interaction in gapped graphene possesses distinctive properties. Our numerical 
results in Fig. [3] showed that although the presence of a gap does not change the ferromagnetic 
nature of the interaction within a chosen sublattice, it introduces an asymmetry between 
these two sublattices. We find that the larger the band gap, the smaller is the ferromagnetic 
energy for one of the two sublattices and the opposite behavior for another sublattice at 
the same time. The antiferromagnetic interaction is also reduced as the gap is increased for 
undoped graphene but not for doped graphene. Moreover, we found the gap modulation 
to the distance dependence of the RKKY interaction for both undoped (1/-R 3 ) and doped 



11 



(mixed 1/R 2 and 1/-R 3 ) graphene, where R is the distance between two magnetic impurities 
on sublattice sites. Additionally, we showed that the doping-induced oscillations in the 
magnetization within a sublattice or between sublattices can be substantially modified by 
an energy gap in both magnitude and phase. Some time ago, Gumbs and Glasser^ showed 
that the RKKY interaction in a metal is modified by the presence of a surface and its long- 
distance behavior is not given by a simple power law as in the bulk. The results in the 
present paper again confirm that the RKKY interaction may be significantly influenced by 
chemical composition or the local environment of the material. 

Acknowledgments 

This research was supported by the contract # FA9453-11-01-0263 of AFRL. DH would 
like to thank the Air Force Office of Scientific Research (AFOSR) for its support. 



1 M. A. Ruderman and C. Kittel, Phys. Rev. 96, 99 (1954). 

2 T. Kasuya, Prog. Theor. Phys. 16, 45 (1956). 

3 K. Yosida, Phys. Rev. 106, 893 (1957). 

4 E. Kogan, Phys. Rev. B 84, 115119 (2011). 

5 B. Uchoa, T. G. Rappoport, and A. H. Castro Neto, Phys. Rev. Lett. 106, 016801 (2011). 

6 M. Sherafati and S. Satpathy, Phys. Rev. B 83, 165425 (2011). 

7 A. M. Black-Schaffer, Phys. Rev. B 81, 205416 (2010). 

8 L. Brey, H. A. Fertig, and S. Das Sarma, Phys. Rev. Lett. 99, 116802 (2007). 

9 S. Saremi, Phys. Rev. B 76, 184430 (2007). 

10 V. K. Dugaev, V. I. Litvinov, and J. Barnas, Phys. Rev. B 74, 224438 (2006). 

11 M. A. H. Vozmediano, M. P. Lopez-Sancho, T. Stauber, and F. Guinea, Phys. Rev. B 72, 
155121 (2005). 

12 P. K. Pyatkovskiy, J. Phys.: Condens. Matt. 21, 025506 (2009). 

13 S. Y. Zhou, G-H. Gweon, A. V. Fedorov, P. N. First, W. A. de Heer, D-H. Lee, F. Guinea, A. 
H. Castro Neto, and A. Lanzara, Nat. Mater. 6, 770 (2007). 

14 G. Li, A. Luican and E. Y. Andrei, 2008 |arXiv:0803.4016| 



12 



O. V. Kibis, Phys. Rev. B 81, 165433 (2010). 

O. Roslyak, G. Gumbs, S. Mukamel, J. Chem. Phys. 136, 194106 (2012). 

A. Iurov, G. Gumbs, O. Roslyak, and D. H. Huang, J. Phys.: Condens. Matt. 24, 015303 (2012). 
D. Huertas-Hernando, F. Guinea, and A. Brataas, Phys. Rev. B 74 155426 (2006). 
Y. Yao, F. Ye, X-L. Qi, S-C. Zhang, and Z. Fang, Phys. Rev. B 75, 041401 (2007). 
G. Gumbs and M. L. Glasser, Phys. Rev. B 33, 6739 (1986). 
Bessel functions are special cases of the general Meijer functions 



13 




FIG. 1: (Color online) Exchange interaction energies Vi, V2 and V3 (left panel), as well as the 
changes of them AVi , AV2 and AV3 (right panel) , in undoped graphene as functions of the half- 
gap Eg for fixed Raa = Rbb = Rab = 0.5. Here, we measure energy in units of 70 and distance 
in units of R^ u = | r M — r' u \ ja for fi, v = A, B. 
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FIG. 2: (Color online) Exchange interaction energies R? AA V\ (top panel), R BB V2 (middle panel) 
and R\bP% (bottom panel) on the left column, as well as the changes of them R\ A AVi (top panel), 
R BB AV2 (middle panel) and R? AB AVz (bottom panel) on the right column, in undoped graphene 
as functions of the distance between two magnetic impurities. Curves 1, 2, 3 correspond to gap 
energy parameter E g = 0.1, 0.3, 0.5, respectively. Here, we measure energy in units of 70 and 
distance in units of R„ 



|r M - r'J/a for /z, v = A, B. 
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FIG. 3: (Color online) Exchange interaction energies Vi, V2 and V3 (left panel), as well as the 
changes of them APi , AP2 and AP3 (right panel) , in doped graphene as functions of the half-gap 
Eg for fixed Raa = Rbb = Rab = 1 and Ep — E g = 1. Here, we measure energy in units of 70 
and distance in units of = Ir^ — r' v \/a for fj,, u = A, B. 
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FIG. 4: (Color online) Exchange interaction energies V\ (top panel), Vi (middle panel) and Vz 
(bottom panel) on the left column, as well as the changes of them AVi (top panel), A 7^2 (middle 
panel) and AV3 (bottom panel) on the right column, in doped graphene as functions of the Fermi 
energy Ep — E g for fixed Raa = Rbb = Rab = 1-0. Curves 1, 2, 3 correspond to the gap energy 
parameter E g = 0.1, 0.3, 0.5 respectively. Here, we measure energy in units of 70 and distance in 
units of R^y = I — r' v \ /a for fi, v = A, B. 
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FIG. 5: (Color online) Exchange interaction energies R AA Vi ( to P)> r2 bb v 2 (middle) and R 2 AB V3 
(bottom) on the left column, as well as the changes of them R 2 AA /S.V\ (top), R BB AV2 (middle) 
and R? AB AVz (bottom) on the right column, in doped graphene at Ep — E g = 1 as functions of 
the distance between two magnetic impurities. Curves 1, 2, 3 correspond to gap energy parameter 
Eg = 0.1, 0.3, 0.5, respectively. Here, we measure energy in units of 70 and distance in units of 
Rfiu = - r'J/a for fj,, v = A, B. 
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